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ABSTRACT 

We studied the decoupling of photons from uhra-relativistic sphericahy symmetric outflows expand- 
ing with constant velocity by means of Monte-Carlo (MC) simulation. For outflows with finite width 
we confirm the existence of two regimes: photon thick and photon thin introduced recently by Ruffini, 
Siutsou, Vereshchagin (2011), hereafter RSV. The probability density function of photon last scatter- 
ing is shown to be very different in these two cases. We also obtained spectra as well as light curves. 
In photon thick case, the time integrated spectrum is much broader than the Planck function and its 
shape is well described by the fuzzy photosphere approximation introduced by RSV. In the photon 
thin case we confirm the crucial role of photon diffusion, hence the probability density of decoupling 
has a maximum near the diffusion radius, well below the photosphere. Its spectrum has Band shape. 
It is produced when the outflow is optically thick and its peak is formed at diffusion radius. 
Subject headings: Gamma Ray: bursts — methods: numerical — radiation mechanisms: thermal 
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1. INTRODUCTION 

Photospheric emission is considered among the lead- 
ing mechanisms generating gamma rays from relativis- 
tically e xpanding plasma in cosmic Gamm a Ray Bursts 
(GRBs) (IGoodmanlll986D. In both firebaU (|Piran II 19991) 
and fireshell (|Rufltini I I2009D models photospheric emis- 
sion is naturally produced when expanding initially op- 
tically thick plasma becomes transparent for photons. 

While within the fireball model the basic approxima- 
tion involved is that of (infinitely long) steady wind, in 
the fireshell model th e geom etry of the outflow is a thin 
shell It is shown by IRS VI that in both wind and shell 
models when outflows with finite duration are considered 
some important aspects of the optical depth behaviour 
are overlooked. There a new physically motivated clas- 
sification was introduced: photon thick and photon thin 
outflows with respectively Rph < Rt and Rph > Rt with 
Rph being the photospheric radius, Rt = 2T^l the radius 
of transition, F is the Lorentz factor and / the labora- 
tory width of the outflow. Light curves and spectra from 
finite outfiows were obtained using two different approx- 
imations: sharp and fuzzy photosphere. These approxi- 
mations were derived there from the transport equation 
for radiation. 

In this work we examine the validity of these approxi- 
mations using MC simulations of photon scattering. We 
compute the probability of last scattering in finite rela- 
tivistic outflows. Then we obtain light curves and spectra 
from both photon thick and photon thin outflows, and 
compare them with the results obtained before. 

2. THE MODEL 



Following IRS VI we consider a spherical wind of finite 
width I with laboratory number density of electrons be- 
ing zero everywhere except the region R{t) < r < R{t) + l 
where it is 
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with Rq and no being respectively the radius and the 
density at the base of the outfiow. This outfiow expands 
with constant Lorentz factor F, so that R{t) — Rq + 

I3ct, where /? = Vl — F"^. Accelerating outflows are 
not considered here. Notice that when there is strong 
relative difference of Lorentz factors within the outflow 
with AF ~ F a radi al spr eading phenomenon disc ussed in 
iPiran et al~l (|1993D and iMeszaros et al. I (|1993D occurs. 
The outflow width starts to increase linearly with R{t) 
at the radius 



Rs = F"; 



F 

AT ■ 



(2) 



It is clear from this formula that such phenomenon occurs 
at radii much larger than T^l when AF ^ F. In what 
follows we consider such case only. Our results are valid 
for arbitrary constant I. The expression for the o ptical 
depth in the model defined by Eq. (1) is found in IRSVI 
and it reads: 



T(r, 9, t) ^ tqRq 
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and r is the position of photon emission, 9 is the angle 
between momentum of photon and the radius vector, a 
is the Thomson cross section, t is the time the photon 
remains within the outflow, Eq is total energy release 
in the GRB, nip is the proton mass, and c the speed of 
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light. Time t is found via the equations of motion of the 
outflow and that of the photon. 

The optical depth in finite relativistic outflows has two 
different asymptotic, depending on initial conditions: in 
photon thick case t is almost constant within the outflow; 
instead in the photon thin case it is linearly increasing 
with dept h from the outer boundary at R{t) + I. 

In IRS VI diffusion of photons within the outflow is found 
to play a crucial role in photon thin case: the radius at 
which photons effectively diffuse out of the outflow is 
smaller than the photospheric radius Rph- The latter is 
defined by equating expression ([3]) to unity for emission 
with = and t corresponding to the inner boundary. 
Hence, when the photon thin outflow reaches the photo- 
spheric radius there are few photons left in it. 

3. MONTE CARLO SIMULATION 

We have used a spherically symmetric MC simulations 
of photon scattering inside the outfiow described by Eq. 
([1}. For more detail on interactions in relativistic ex- 
panding plasma (see Begue, Siutsou and Vereshchagin, in 
preparation). In this simulation each photon is followed 
as it interacts with electrons until it decouples from the 
outflow. Photons are injected in the outflow when the op- 
tical depth of the inner boundary is given in Tab. [T] It 
is also assumed that their initial distribution is isotropic 
and thermal. The radial position of each photon inside 
the outfiow is chosen randomly. 

The code consists mainly of a loop computing each 
scattering. We proceed in two step s. First an infinite 
and steady wind (already treated by Belo borodov I ()2011[ ) 
with / — >■ oo and i? — >■ 0) is considered. For a given po- 
sition characterized by r and 9 of photon in the outflow, 
we compute a maximal value for the optical depth Tmax 
using ([3]) with t — > oo. The probability for the photon 
to decouple the outfiow is exp{—Tmax)- Then a ran- 
dom number X £ [0, 1] is chosen. On the one hand, if 
X < exp (—Tmax) the photon is considered decoupled. 
Afterwards photon remains in the outflow but does not 
scatter. This case corresponds to decoupling in photon 
thick case, for which the presence of boundaries is not 
essential. On the other hand, if X > exp {—Tmax) the 
position 7's of next scattering is computed from the opti- 
cal depth 
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where 9s is the angle between photon momentum and the 
radius vector at the position of scattering. The new radial 
position Ts such that exp [— r(r, Ts, 0)] = X is found by 
iterations. 

The second step is to take into account the flnite width 
of the outflow. To do so, Ts is compared with the radii 
of the inner and outer boundaries of the outflow. If < 
R{ts) or Vs > R{ts) + l the scattering does not take place 
and the photon is considered decoupled. Such decoupling 
occurring at the boundaries corresponds to the photon 
thin case. In the opposite case R{ts) < < R{ts) + l the 
scattering is assumed to occur and the loop is repeated 
until the photon decouples. 

We consider two models of scattering: coherent 
isotropic scattering and Compton scattering. The for- 
mer is treated for comparison with results in the litera- 
ture such as e.g. iBeloborodovl (|2011[) . Such model is also 



interesting per se since it preserves Planck spectrum and 
it traces only geometrical effects. Instead when Comp- 
ton scattering is considered, the equilibrium spectrum is 
the Wien one because stimulated emission is not taken 
into account by MC simulation. For initially large opti- 
cal depth the spectrum indeed firstly acquires the Wien 
shape, and only at the photosphere it changes the shape 
again. 

The coherent scattering is computed by Lorentz trans- 
formation to the reference frame comoving with the out- 
flow. In addition, when Compton scattering is consid- 
ered, another Lorentz transformation to the rest frame 
of the electron is performed. The electron is chosen ran- 
domly from the Boltzmann distribution at a given co- 
moving temperature defined as 



T = TqT 
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where a = 4c r<;B/c, asB is the Stefan-Bol tzmann con- 
stant, see e.g. iRuffini, Vereshchagin I (|20ll) . 

Various tests of the code are performed. Before deal- 
in g with finite outflow s we also reproduced the results 
of iBeloborodov I (|2011[) . We present below the results 
for two specific cases. The set of parameters for these 
simulations are given in Tab. [TJ The left (right) column 
corresponds to photon thick (thin) case. 



TABLE 1 
Parameters for the simulations. 





Photon thick 


Photon thin 


-Ro (cm) 


IQS 


10» 


I (cm) 


4.5 * 10* 


10» 


TO 


2.3 * 10^" 


1.2 * 10" 


Eo (erg) 


1.5 * 10^^ 


lOS-i 


r 


500 


300 


Rph (cm) 


4.6 * 10^^ 


3.3 * 10^"* 


Rd (cm) 


3.5 * 10" 


1.0 * 10^* 


n 


10^ 


1.3* 10=* 
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4. RESULTS 

Probability density function of position of last 
scattering 

Consider first the probability density of photon last 
scattering position as a function of depth (top panel of 
Fig. [1]). Photon decoupling from photon thick outflows 
is expected to be local and the presence of boundaries 
should not change this probability substantially, indeed 
the probability distribution function of last scattering is 
found to be almost independent on the depth. On the 
contrary in photon thin outflow there is enough time for 
photons to be t ransp orted at the boundaries by diffusion 
as discussed in IRSVI . As a result the probability density 
is peaked at the boundaries. 

Difference between photon thin and photon thick out- 
flows is also reflected in the probability density of last 
scattering as a function of the radius, shown at the bot- 
tom of Fig. [TJ In the photon thin outflow most photons 
escape from the outfiow well before the photospheric ra- 
dius, namely near diffusion radius 



Rd = {toT^RqI' 



,1/3 



(7) 



MC simulation of photospheric emission 
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Fig. 1. — Upper panel: Probability density of last scattering as 
a function of normalized depth for photon thick and photon thin 
outflows (in the latter case decreased by a factor of 10 for a better 
presentation). 

Lower panel: probability density function for the position of last 
scattering in the following cases: infinite and steady wind (1), pho- 
ton thick case of Tab. [T] (2), photon thick case of Tab. [l] with a 
different Lorentz factor F = 300 (3), and photon thin outflow (4). 
The vertical line (a) represents the diffusion radius in photon thin 
case, line (b) represents the photospheric radius while lines (c) and 
(d) show the transition radius Rt for curves (2) and (3), respec- 
tively. Both coherent and Compton scattering models give the 
same results. 



Instead for photon thick outflows the probability den- 
sity function of last scattering is found to be close to the 
one of infinite and steady wind, found by iBeloborodov I 
(|2011| ) . The finite extension of the outflow results in the 
exponential cut-off for the probability density function 
at radii larger than Rt- 

Let us stress how different the probability densities for 
photon thick and photon thin cases are. Not only the 
positions of the maximum, but also the shapes are dif- 
ferent. In photon thin case at small radii the number of 
photons that diffuse out is determined by the change in 
the diffusion coefficient, i.e. the electron density which 
follows a power law. In expanding plasma while the den- 
sity decreases, the mean free path for photons increases 
which makes the diffusion faster. At large radii almost 
all photons have already diffused out, and the probability 
has an exponential cut off. As for the photon thick case, 
the probability density dependence on the radius is oppo- 
site. For small radii the probability for last scattering is 
exponentially suppressed because of large optical depth, 
w hile at l arge r adii it follows a power law, as discussed 
in iPe'er I ()2008D . At larger radii we find an exponential 
decrease due to the finiteness of the outflow, see Fig. [T] 
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Fig. 2. — Average number of scatterings as function of initial 
optical depth r^. Results of MC simulations are shown with points, 
while analytic results in different asymptotic from Appendix are 
shown by the corresponding lines. In the optically thin regime the 
number of scatterings is proportional to r^. For 1 < < r(i? /j) it 
is proportional to t^. For t(R0) < Ti < T{Rt) it is proportional 
1/2 

to ' . For even larger Ti corresponding to the photon thick 
asymptotic the number of scatterings is again proportional to Ti. 
In these simulations we used F = 100, I = 10* cm and tq = 10^^. 

We also found the difference between photon thick and 
photon thin cases when computed the average number 
of scatterings as function of initial optical depth, shown 
in Fig. [2] Derivation of analytic results is reported in 
Appendix. 

4.2. Spectra and light-curves from photon thick outflow 

The average comovin g ener g y of phot ons for steady 
wind was considered bv IPe'er I (|2008| ) and IBeloborodov I 
(|2010f l. For finite outflow we show this quantity in Fig. [3] 
and compare it with the r~'^l^ dependence characteristic 
of adiabatic cooling. The difference between Compton 
and coherent scattering comes from the fact that the av- 
erage energy in the optically thick regime is SfcT for Wien 
spectrum and 2.82fcT for Planck spectrum. In agreement 
with previous studies we find that the average energy of 
photons at the photosphere is higher than adiabatic cool- 
ing predicts. It is 1.58 in the case of Compton scattering 
and 1.4 in the case of coherent one. 

In Fig. |3]we compare time integrated spectra of pho- 
ton thick outflow obtained with different models. Each 
of them involves 5 * 10^ photons. In order to compare the 
results of IRS VI with our MC simulations we shifted the 
spectra for sharp, fuzzy approximations and Plank spec- 
trum by the factor 1.58, which takes into account the 
effect discussed above. The results obtained in fuzzy ap- 
proximation are in good agreement with MC simulations. 
The later gives the low energy spectra index a = 0.24 for 
Compton scattering and a — .19 fo r coherent one to be 
compared with a = found in IRS VI . 

We also show in Fig. [S]time resolved spectra at selected 
arrival times which are clearly far from being Planckian. 
T he effec t of sp ectral flattening at late time discovered 
bv IPe'er I (|201ll ) is visible. 

Finally, light curves are presented in Fig. [6l where we 
defined t^h — l/c and ta = t^ — ^Re/c, with r^ and te 
laboratory radius and time of emission. We chose ta = 
for a photon emitted at Re = Rq- The results of MC 
simulation are in qualitative agreement with both ap- 
proximations, slightly favouring the sharp photosphere 
one. 
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Fig. 3. — Average comoving energy of photons in photon thick 
outflow as function of the radius. The averaging is performed for 
photons still coupled with the outflow. 



Fig. 6. — Light curves from photon thick outflow within differ- 
ent approximations: sharp photosphere (1), fuzzy photosphere (2) 
and MC simulation (3) for Compton scattering model. Coherent 
scattering model gives similar results. 




hv/kTLos 

Fig. 4. — Time integrated spectrum of photospheric emission 
from photon thick outflow computed within different approxima- 
tions: Planck spectrum (1), sharp photosphere approximation (2), 
fuzzy photosphere approximation (3), Monte Carlo simulation for 
Compton scattering (4) , Monte Carlo simulation for coherent scat- 
tering (5). Spectra (1), (2) and (3) are shifted in energy by a factor 
1.58 (see text for details). T^os the laboratory temperature of 
the outflow at the photospheric radius. 
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Fig. 5. — Time integrated spectrum from photon thick outflow (1) 
together with time resolved spectra corresponding to the beginning 
of the plateau in the light curve (2), middle of the plateau (3) and 
the tail (4). 



4.3. Spectra and light curves from photon thin outflow 

In photon thin case diffusion has been shown to be 
important. Both hght curve and spectrum are charac- 



terized by the radius of diffusion Ro and the associated 
arrival time tn — Rd/{'2^^c). It should be noted that 
the photospheric emission carries only a fraction of the 
total energy of a GRB. In the photon thin case the du- 
ration of the photospheric emission is relatively short, 
see the formula above, which implies that the luminosity 
of this component can be comparabl e to the luminos- 
ity of the rest of the prompt emission (jBernardini et alH 
I2007f ). Such short and intense photospheric emission pre- 
ceding the main burst has been reported in many GRBs 
(|Rvde and Pe'erl[200l . 

Since photons are still strongly coupled to the outflow, 
local thermodynamic equilibrium is maintained every- 
where but at the very boundary. Indeed, this part of 
the outflow looses photons by random walks. For MC 
simulations of the photon thin outflows we used the co- 
herent scattering model. The results are compared with 
IRS VI where light curve and spectrum are obtained us- 
ing the sharp photosphere approximation, after solving 
radiative diffusion equation. It has been demonstrated 
there that the temperature at the boundary of the out- 
flow follows a different law from the one of adiabatic 
cooling: T{r) oc 7'-i3/24 -^^ith both temperatures nearly 
coinciding at the radius of diffusion. So we expect the 
time integrated spectrum obtained by MC simulations 
has more power at energies above the peak. 

Finally, in both MC simulations and the work of IRS VI 
the photons emitted outside the relativistic beaming cone 
from the photon thin outflows are neglected. This is 
good approximation, since such photons have much lower 
energy due to the Doppler effect. Besides, their arrival 
time is 2V^ longer than the arrival time of diffusion. 

Spectra and light curves are shown respectively in Figs. 
[7]and[8l The simulation involves 10^ photons. We see a 
good agreement between both methods of computation. 
The high energy tail with a power law in the time inte- 
grated spectrum is clearly visible. Our MC simulation 
shows the exponential cut off at the energy correspond- 
ing to the temperature of photons at the radius defined 
by the initial optical depth t; (see Fig. [7]). The power law 
above the peak e nergy should extend up to initial tem- 
perature To, see iRSVl . This also results in the missing 
initial part of the light-curve (see Fig. [5]). 

The power law index is j3 = —3.5 for the model of the 
outflow given by Eq. ([!]). Such power law is steeper 



MC simulation of photospheric emission 
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Fig. 7. — Time integrated spectrum from photon thin outflow ob- 
tained by MC simulation (crosses), and sharp photosphere approx- 
imation (thick curve). Ti^os is here the laboratory temperature at 
Rd- 



•a 



than typical high energy power laws observed in GRBs. 
It remains to be analyzed how (3 changes for different hy- 
drodynamic models of the outflow. Besides, high energy 
GcV emission cannot be explained with this model since 
there exists an exponential cut off corresponding to the 
initial temperature in the source of GRB, typically being 
in MeV region. 



5. CONCLUSIONS 

To summarize, in this work we considered the photo- 
spheric emission from finite relativistic outflow by MC 
simulations of phot on sca ttering. The validity of the as- 
sumptions made in IRSVI is verified and it is found that 
the fuzzy photosphere approximation used in that paper 
is adequate in describing both the light curves and the 
shape of the spectra of photospheric emission from pho- 
ton thick outflows. The sharp photosphere approxima- 
tion produces qualitatively acceptable results for photon 
thick outflows, but underestimates high energy part of 
the spectrum and the decaying part of the light curve. 

As expected, diffusion is found to dominate the pho- 
tospheric emission of photon thin outflows. The proba- 
bility density of last scattering is obtained in this case 
both as function of outflow depth and its radial position. 
Both these functions differ from their photon thick coun- 
terparts, once more pointing out the essential difference 
between these two cases. Good agreement is obtained 
between MC simulation and sharp photosphere approx- 
imation for both light curve and spectrum from photon 
thin outflow. 



Fig. 8. — Light curves from photon thin outflow obtained by MC 
simulation (crosses) and sharp photosphere approximation (thick 
curve) . 
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APPENDIX 

AVERAGE NUMBER OF SCATTERINGS 

We derive here the average number of scatterings expected for photons in both photon thin and photon thick 
outflows. 
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Froni|RS3 we recall the expressions of the optical depth for these two cases: 



2r2 



r 

The average number of scatterings is defined as 



To-i-, r»r2; 



(Al) 



(AT) = / ^ = / a{n,)cdt,, (A2) 



where tc is comoving time, A = 1/ {aTn-c) is the comoving mean free path of photons, Uc is comoving density, ti is the 
initial comoving time, tf is final comoving time when the photon leaves the outflow. The integral (jA2l) should be taken 
along the average photon path. In the case of optically thick medium this path is given by the equation of motion of 
the outflow 

r = r, + Pet, (A3) 

where is the laboratory radial position of the photon at initial laboratory time. 

In the photon thick case the photons stay in the outflow long after decoupling, sotf oo. Then, taking into account 
relations the ric = ni/T, dr = f3cdt, and tc = t/T along the world-line (jA3[) . we obtain 

(N) = a-— = 2r„ (A4) 

where is optical depth of the outflow at r^. This result is in agreement with lVurm et al.l (|2011[ ). 

In photon thin case the photons leave the outflow from its boundaries, and the time interval needed for the photon 
to reach them by random walk can be estimated as fj^ = 1^/ Dc, where Ic — T/ is the comoving radial thickness of the 
outflow, and diffusion coefficient is Dc = (cAc)/3 = c/(3anc). When this time is much less than the characteristic time 
of expansion, equal to ti, that is the case when initial radius r.i is much larger than radius of diffusion Ro, we have 

{N) 0, irl (A5) 

In the opposite case when initial radius is much smaller than the radius of diffusion Rd, we have 



W^^\jnr,^. (A6) 



